William I. Fine Theoretical Physics Institute 
University of Minnesota 



FTPI-MINN-13/09 
UMN-TH-3141/13 
March 2013 



Mixing of partial waves near B*B* threshold in e + e 

annihilation 

Xin Li a and M.B. Voloshin a fe c 

'School of Physics and Astronomy, University of Minnesota, Minneapolis, MN 55455, USA 
b William I. Fine Theoretical Physics Institute, University of Minnesota, 
Minneapolis, MN 55455, USA 
institute of Theoretical and Experimental Physics, Moscow, 117218, Russia 



Abstract 

We consider the production of B*B* meson pairs in e + e~ annihilation near the 
threshold. The rescattering due to pion exchange between the mesons results in a 
mixing between three partial wave amplitudes: two P-wave with the total spin of 
the meson pair S = and S = 2 and an P-wave amplitude. The mixing due to 
pion exchange with a low momentum transfer is calculable up to cm. energy E ~ 
15 -T- 20MeV above the threshold. We find that the P — F mixing is numerically quite 
small in this energy range, while the mixing of the two P-wave amplitudes is rapidly 
changing with energy and can reach of order one at such low energies. 



The strong interaction between hadrons containing heavy c or b quarks as well as the light 

ones results in a peculiar dynamics at energies where these hadrons move relatively slow with 

respect to each other. In particular, the cross section for production of heavy meson pairs 

in e + e~ annihilation displays an intricate behavior near the thresholds for new meson pair 

channels. This behavior, observed for the onset of the production of charmed mesons pQas 

well as of the heavier B and B* meson pairs [21 [3], is not yet known in detail and further 

studies may uncover new hadronic structures at these thresholds and provide an insight into 

the dynamics of systems containing both heavy and light quark degrees of freedom. One 

important yet unknown feature of the production of pairs of heavy vector mesons is the 

spin-orbital structure of the production amplitudes. Namely, the overall quantum numbers 
JPO = i f 

a pair of e.g. B*B* mesons produced in the e + e annihilation allow for three 
different combinations of the orbital momentum and the total spin S of the pair: P wave 
with S = as well as a P or F wave with S = 2. It is generally expected that the F wave is 
kinematically suppressed near the threshold, which still leaves unknown the composition of 
the P-wave production amplitude in terms of the S — and 5 = 2 components. The latter 
composition, clearly measurable from angular distributions [I], can be quite nontrivial [5] 
and in fact rapidly varying function of the cm. energy in the near threshold region. In 
either case, the actual composition of the amplitude is very likely to be much different from 
a naive expectation [SJ El |4] of dominance (by a factor 20 in the cross section) of the S = 2 
wave, following from the heavy quark spin symmetry, which symmetry is naturally broken 
in the threshold region 

It is quite clear that the forces between the heavy mesons depending on the spins of light 
quarks result in a mixing of the partial waves, as recently discussed [8] for the S — D mixing 
in the J p = 1 + channel, and for heavy mesons the effect of these forces is enhanced by the 
factor of the meson mass M. There is likely a part of this interaction at short distances 
determined by Aqcd that is currently impossible to analyze in a model independent way. 
However the interaction also contains a long-distance part due to the pion exchange (see 
e.g. in Ref. [9]), determined by the strength g of pion interaction with the heavy mesons, 
known [TO] from the decay rate D* — > Dir. The effects of the pion exchange can be separated 
from those of the short-distance interaction in the range of cm. energy above the threshold 
E = p 2 /M where the cm. momentum p is small compared to Aqcd, P 2 Aq CD . Indeed, 
in the J PC = 1 channel two types of effects are possible: the P — F wave mixing and the 
mixing of the S = and S = 2 P waves. The short-distance contribution to the P — F wave 
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mixing is determined by the parameter Mp 2 /Aq CD , while the pion contribution behaves [8] 
as g 2 M/Aqcd at p 2 ^> /x 2 and as g 2 M p 2 /(AqcdH 2 ) when p is comparable with or smaller 
than the pion mass \l. In what follows we calculate the effects arising from the pion exchange 
in the partial wave mixing for pairs of heavy vector B* mesons^, so that in our numerical 
estimates we use M = 5325 MeV. We find that the P — F mixing, although parametrically 
enhanced at low p, is still of a moderate value and reaches only about 0.1 in the amplitude 
at the upper end of the applicability of our approach, p pa 300 MeV. Thus we confirm the 
existing expectation that the presence of the F wave in production of heavy vector meson 
pairs in e + e~ annihilation is likely insignificant at energies slightly above the threshold. 

The mixing of the S = and S = 2 channels due to the short- distance interaction is 
generally of order one at any energy near the threshold. However, the energy scale for a 
variation of this part is set by Aq CD /M, so that no significant change is expected as long 
as p 2 <C Aq CD . In particular, the absorptive part of the S = and S = 2 mixing amplitude 
is proportional to M p 3 /Aq CD due to the P wave phase space. On the other hand, the pion 
exchange contribution experiences a significant variation on a smaller energy scale. Namely 
we shall argue that the absorptive part of the S = and S = 2 mixing amplitude behaves as 
g 2 M p/ Aq CD at fi 2 <C p 2 <C Aq CD and numerically changes from zero at the threshold to a 
factor of order one at p ~ 200 -f- 300 MeV. Thus the expected effect of the pion exchange in 
the latter mixing is a rapid variation above the threshold in the range of excitation energy 
up to 15 ^20 MeV. 

It can be also noticed that the re-scattering between the channels with two vector mesons 
and those with one or two pseudoscalar mesons, e.g. BB — > B*B* and B*B — > B*B*, which 
generally also contributes to the mixing of partial waves of the vector mesons, should be 
considered as a short-distance effect, whether it proceeds through the pion exchange or 
through other forces. Indeed, the momentum transfer in these processes is of the order of 
q ~ vMA with A being the mass difference between the vector and pseudoscalar mesons. 
Since in the heavy quark limit the parametric behavior is A ~ Aq CD /M, one finds that in 
the cross-channel re-scattering q ~ Aqcd, and at such momentum transfer the pion exchange 
is indistinguishable from other short-range forces. For this reason in our calculations of the 
pion exchange we consider only the diagonal re-scattering B*B* — > B*B* as shown in Fig. 1 

The diagonal interaction of pions with the isotopic doublet of heavy vector mesons V = 

1 A similar calculation is applicable also to the D*D* meson pairs. However the significance of the 
discussed effects at a fixed cm. momentum is scaled down by the lighter mass of the charmed mesons. 
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Figure 1: Diagonal re-scattering through pion exchange of B*B* meson pairs produced in 
e + e~ annihilation 



(V + , V°) can be written as 



H int = tfe m (V J U a V k )d l n a 



where r a are the isotopic Pauli matrices, f n ~ 132 MeV is the pion decay constant, and g is 
a dimensionless coupling constant determined from the D* — > Dtt decay as g 2 ~ 0.18. It is 
assumed in Eq. (DO) as well as throughout this paper that the nonrelativistic normalization is 
used for the wave functions of the heavy mesons. 

The amplitude for the production of a B*B* pair in e + e~ annihilation near the threshold 
can generally be written in term of three partial wave amplitudes: 



1 1 / 2 

A(e + e~ B*B*) = A (p 2 ) j k p k ■ - afti + j k ■ -j= [c^bj + - -dij a t bi ) x 

lA 2 (p 2 )6 kiPj + ^=A F (p 2 ) 



-^PiPjPk - \ ipk^ij +PjS ik +PiS kj ) 
p o 



(2) 



where p is the cm. momentum of one of the mesons, a and b are the polarization amplitudes 
for the meson and anti-meson and j denotes the polarization amplitude of the virtual photon. 
The amplitudes Aq and A 2 are the S = and S = 2 P-wave amplitudes and Ap is standing 
for the F-wave one. The relative normalization of the amplitudes in Eq. ([2]) is chosen in such 
a way that the production cross section is proportional to p 3 (|^4o| 2 + I A 2 \ 2 + \Ap\ 2 ). One can 
also notice that under this normalization the expansion of Ap at small momentum p starts 
with p 2 . 
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In what follows we treat the mixing of partial waves induced by the pion as a small 
effect and we calculate it in the first order of perturbation theory for which we use the 
nonrelativistic (in heavy mesons) formalism. Proceeding in this way and considering the 
projection on the F-wave we find the following expression for the amplitude A F generated 
after re-scattering through the pion exchange by the P-wave amplitudes A and A 2 . 



A F (p 2 ) 
d 3 q 



I 



(2nf 



£1 

f 2 P 2 

/ 30 



1 PiPjPk - 7 (pk$ij +Pj8ik +Pihj) 



p* " " 5 
A (q 2 ) + A 2 (q 2 ) 



x 



_^ qi(qj-Pj) {qk-Pk) 

■ p 2 — ie (q — p) 2 + fi 2 



(3) 



One can notice that this expression includes the isotopic factor of 3, which corresponds to 
the pion exchange interaction in the isoscalar state of the B*B* pairs produced in e + e~ 
annihilation. The presence of the F wave projector in Eq.([3]) implies that only the part of 
the integral proportional to PiPjPk contributes to the P — F mixing. Namely, if one writes 
the general expression allowed by the symmetry for the integral 



d?q 



A (q 2 ) + A 2 (q 2 



M 



Qi (Qj ~Pj) (Qk-Pk) 



q 2 — p 2 



ie 



p) 2 + fi 2 



Ci(p 2 ) Pi $jk + C 2 {p 2 ) (jpj 6 ik + p k Sij) + C 3 (p 2 ) PiPjPk , 
only the structure proportional to C3 contributes to the expression 

2g 2 



(4) 



A F (p 2 



5f 2 



p 2 C 3 (p 2 



for the amplitude A F : 

(5) 



Clearly, if the amplitudes Aq and A 2 are smooth functions as q 2 varies on the scale of p 2 
or fi 2 the integral for C3 converges and is determined by the range of q such that q 2 ~ p 2 , or 
q 2 ~ fi 2 if p 2 < fi 2 . In order to estimate the numerical significance of the P — F mixing we 
approximate the amplitudes Aq and A 2 by constants, in which case the integral is calculated 
analytically and the result reads as 



A F (p 2 ) = r FP (p) 
with the mixing function r F p(p) given by 



A + A 2 



(6) 



r FP (p) 



g 2 M/i (5 + 18t + 8t 2 



+i 



207T f 2 

15 + 24t-4t 2 
24 1 3 / 2 



16i 5/2 

5+ 18t + 8t 2 



32t 5 /2 



arctan(2\/t) 
log(l + 4t) 



15 + 34t 
24t 2 



(7) 
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where t = p 2 /pi 2 . 

The plot for the function rpp(p) is shown in Fig. 2. One can see that the discussed 
P — F wave mixing is quite small. Indeed, at the upper end of the applicability range of 
our calculation, at p « 300 MeV, the mixing function is still less than 0.05, corresponding 
to only of order 0.1 mixing with the S = P-wave amplitude, and smaller for the S = 2 
P-wave amplitude. 
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Figure 2: The function rpp describing the mixing of the P and F waves: the absolute value 
(solid), and the real (dotdashed) and imaginary (dashed) parts. 



The available data [21 E] do not indicate a presence of a resonance at the B*B* threshold 
in e + e~ annihilation. However the data are not yet conclusive and a threshold resonance may 
show up in future studies. In this case our approximation of smooth (constant) amplitudes 
Aq and A2 would generally be not applicable literally. For this reason we have verified that 
using in Eq.([3]) these amplitudes with a Breit-Wigner shape instead of constants does not 
qualitatively change the conclusion that the P — F mixing is small and remains at the level 
of 0.1 or less for a broad range of the resonance parameters. 

We shall argue that the effect of the pion exchange is significantly larger numerically for 
the mixing of the two P-wave amplitudes Aq and A^. Proceeding to a calculation of this 
effect we notice that the mixing Aq — > A2 is the same as A2 — > Aq by reversibility, so that 
it is sufficient to consider the mixing only 'in one direction', e.g. the S = 2 amplitude A 2 
generated by the S = production amplitude Aq. Using Eq.([T]) for the heavy meson - pion 
interaction and our definition in Eq.([2]) of the production amplitudes, we find that for a pair 
initially produced by the amplitude Aq an S = 2 state is generated by rescattering with the 



5 



amplitude 



A n2 = jk 
V5g 2 



1 



20 
d 3 q 



ciibj + djbi - ^5ij aibi ) x 



M 



Qk (qi-Pi) (qj-Pj) 



f 2 J (27r) 3 q 2 — p 2 — ie (q — p) 2 + fi 2 

If one writes the general expression for the three-index integral in terms of partial waves, 



d 3 q 



Ao(q 2 



M 



qk {qi ~ Pi) (qj - Pj) 



q z — p z 



ie 



(q-p) 2 + /i 2 



D 1 (p 2 )p k 5 i 



-D 2 (p 2 ) ( Pj S ik + Pi 5 jk ) + D 3 (p 2 



\piPjPk - ~ (Pk^ij +Pj$ik +Pi8 k j) 
p o 



(9) 



it can be readily seen that the S = 2 projector in Eq.flH]) selects only the part proportional 
to the invariant function D 2 . As a result the expression for the generated S = 2 amplitude 
SA 2 can be written as 



6A 2 (p 2 



d 3 q 

V$f 2 J 



A (q 2 



M 



q* — p" — ie p 



1 3(p-q — p 2 
1 



p-q)- ip-q){q-p)' 



(q — p) 2 + fi 2 



(10) 

Unlike the integral for the P — F mixing in Eq.flH]) this expression does not converge for a 
constant A (q 2 ) and thus is not determined by the intermediate momentum q 2 ~ p 2 Aq CD 
if the amplitude A (q 2 ) varies at the scale of Aq CD . Thus the pion exchange at small 
momentum transfer does not dominate the discussed mixing of the P waves, and one should 
take into account other interactions at short distances. However the significance of the mixing 
generated by the pion exchange at longer distances can still be evaluated from Eq. flTOl) by 
considering the absorptive part of the mixing determined by q 2 = p 2 . The calculation of 
the absorptive part of the mixing amplitude is done by replacing in Eq. flTO]) the propagator 
(q 2 — p 2 — ie)^ 1 with rr 5(q 2 — p 2 ), and one finds 



SA 2 (p 2 



abs 



r 20 (p) A (p 2 



111 



with the absorptive mixing function r 20 given by 

- 2 Mp 



r 20 (p) 



9 



16\/5 7r f 2 



Ap 4 



log 1 + 



Ap 2 
1* 



121 



The plot of the function r 2 Q is shown in Fig. 3. One can see that this function changes 
between zero at the threshold to rather large values of about 0.75 at the upper end of the 
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applicability of our calculation. This significant and rapid variation of the mixing in fact 
justifies a consideration of the absorptive part alone, since the dispersive part and any other 
effects arising from short distances are expected to exhibit a variation on the momentum 
scale of order Kqcd-, which scale is parametrically larger than the range of the plot in Fig. 3. 
Thus any possible cancellation in the absorptive part due to the short- distance processes 
cannot take place in the entire range of momenta, at which our approach is applicable. 




Figure 3: The function r 2 o describing the absorptive part of the mixing between the S = 
and S = 2 P- waves. 

In summary. We have considered the effect of the pion exchange on the mixing between 
three partial waves of B*B* mesons produced in e + e~ annihilation at energy near the thresh- 
old. The pion exchange is calculable and dominates the mixing of P and F waves as long as 
the cm. momentum p of the mesons is small as compared to Aqcd, which restricts the range 
of the excitation energy of the meson pair to at most E m 15 -r- 20MeV. We find that the 
P — F mixing is rather small and should not exceed approximately 0.1 in the amplitude. The 
mixing effect is however significantly larger numerically for the mixing of the two P-wave 
states, corresponding to the total spin of the meson pair S = and S = 2. Only for the 
absorptive part of this mixing the dominance of the pion exchange can be ensured, while the 
full effect generally depends on the unknown interaction at short distances determined by 
Aqcd- We find that in the latter case the absorptive part of the mixing rapidly changes with 
energy from zero at the threshold to about 0.75 at the upper end of the range where our 
approach is applicable. We thus conclude that the partial wave composition of the produced 
pairs of vector B* mesons should exhibit a nontrivial behavior near the threshold in e + e~ 
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annihilation, which can be studied experimentally. 
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